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Abstract. We have investigated the effects of plasma instabilities driven by rapid e ± pair cascades, which arise in the 
environment of GRB sources as a result of back-scattering of a seed fraction of the original spectrum. The injection of 
e ± pairs induces strong streaming motions in the ambient medium. One therefore expects the pair-enriched medium ahead 
of the forward shock to be strongly sheared on length scales comparable to the radiation front thickness. Using three- 
dimensional particle-in-cell simulations, we show that plasma instabilities driven by these streaming e ± pairs are responsible 
for the excitation of near-equipartition, turbulent magnetic fields. Our results reveal the importance of the electromagnetic 
filamentation instability in ensuring an effective coupling between e ± pairs and ions, and may help explain the origin of large 
upstream fields in GRB shocks. 
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RPIC SIMULATIONS 

Particle-in-cell (PIC) simulations can shed light on the physical mechanism of particle acceleration that occurs in the 
complicated dynamics within relativistic shocks. Recent PIC simulations using injected relativistic electron-ion and 
electron-positron jets show that acceleration occurs within the downstream jet [1, 2, 3]. In general, these simulations 
have confirmed that relativistic jets excite the Weibel instability. The Weibel instability generates current filaments and 
associated magnetic fields [4], and accelerates electrons [1, 2, 3]. 

Three simulations were performed using an 85 x 85 x 640 grid with a total of 380 million particles 
(27 particles/cell/species for the ambient plasma) and an electron skin depth, Ac e = c/(Op e = 9.6A, where 
U)pe = (47 te 2 n e /m e ) 1 / 2 is the electron plasma frequency and A is the grid size. In all simulations jets are in- 
jected at z = 25 A in the positive z direction. In all simulations radiating boundary conditions were used on the planes 
at z = 0,Zmax* Periodic boundary conditions were used on all other boundaries. The ambient and jet electron-positron 
plasma has mass ratio m e /m v = m e -/m e + = 1. The electron thermal velocity in the ambient plasma is v t h = 0.1c 
where c is the speed of light. The mass ration on ions in the ambient plasma is mj/m e = 20 

In this report we simulate three different initial jet electrons (positrons) distributions. The first case (A) is as same 
as one of the simulations in previous papers [2, 3] (yVj jj = 12.57c). The second (B) and third cases (C) are a special 
case and corresponds to cold jet electrons and positrons created by photon annihilation with 2c < yV j| < 30c [5]. For 
all cases the temperature of jet particles are very cold (0.01c in the rest frame). In this letter we concentrate on the 
simulation with the broad 2 < ypj| < 30 Lorentz factor distribution. The first (A) and second (B) cases are injected 
into electron-ion ambient plasma, the third (C) case is injected in ambient pair plasma. The average jet velocity for the 
broad pair loading is nearly yVjj = 12.5c. 

Figure 1 shows snapshots for three cases (top: A, middle: B, and bottom: C). The Weibel instability in the electron- 
positron jet injected at z/ A — 25 is illustrated in 2D images in the x — z plane (y = 43 A) in the center of the jet at 
t = 59.8/co^e. In the left column the colors indicate the z-component of the current density (7 Z ) with 7 Zj / x indicated by 
the arrows. Behind the jet front (z/ A = 600) the Weibel instability perturbs the electron density, leading to nonuniform 
currents and transverse magnetic fields [2, 3]. The magnetic field energy density averaged in the x — y plane along 
the z-direction. The total energy density (0.5 * (B\ -F By + B ^)) is plotted by the red curves, the perpendicular energy 
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FIGURE 1. The snap shots of simulation with at / = 59.8/cty e for three cases. In the left column the colors indicate the Z- 
component of the current density (7 Z ) with 7 Z ,/ X indicated by the arrows. In the right column the total energy density (0.5 * (fi\ + 
By + B\)) is plotted by the red curves, the perpendicular energy density (0.5 * (B\ -f #y)), by green, and the parallel energy density 
(0.5 * B\), by blue curves. 


density (0.5 * (B\ + B y)), by green, and the parallel energy density (0.5 * Z^), by blue curves. The Weibel instability 
in the electron-positron jet injected at z/A = 25 is illustrated in 2D images in the x — z plane ( y = 43A) in the center 
of the jet at t = 59.8/ 0)p e . In the left column the colors indicate the ^-component of the current density ( J z ) with 7 Z ,7 X 
indicated by the arrows. Behind the jet front (z/ A = 600) the Weibel instability perturbs the electron density, leading 
to nonuniform currents and transverse magnetic fields. 

The comparison among three cases shows that the electron-ion ambient plasma enhances the generation of magnetic 
field due to the excitation of ion filamentation instability. The pair jet with broad distribution sustains the magnetic 
field strength due to the longer linear stage with jet particles with larger Lorenz factor. Further discussions will be 
reported separately [6]. 
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